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QCD with 2 massless colour-sextet quarks is studied as a model of Walking Technicolor We 
simulate lattice QCD with 2 light color-sextet staggered quarks at finite temperature, and use the 
dependence of the coupling at the chiral transition on the temporal extent, A^,, of the lattice in 
lattice units to study the running of the bare lattice coupling with lattice spacing. Our goal is 
to determine whether this theory is QCD-like and 'walks', or if it is conformal. If it is QCD- 
like, the coupling at the chiral transition should tend to zero ss, Nt ^ °° m a. manner controlled 
by asymptotic freedom, i.e. by the perturbative j3 -function. On the other hand, if this theory is 
conformal, this coupling will approach a non-zero limit in the Nt ^ °° limit. We are extending 
our simulations on an A', = 8 lattice to determine the position of the chiral transition with greater 
accuracy, and are performing simulations on an A', = 12 lattice. 
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1. Introduction 

We are interested in extensions of the Standard Model with strongly-interacting Higgs sectors. 
The most promising theories of this type are Technicolor models [jl], ^. Technicolor theories are 
QCD-like theories where the techni-pions play the role of the Higgs field, giving masses to the W 
and Z. Walking Technicolor theories, where the fermion content is such that the coupling constant 
evolves very slowly - 'walks' - and their extensions can avoid phenomenological problems, which 
otherwise afflict Technicolor models [^, ^, ||, ^. 

QCD with 1^ < A'/ < 3^ flavours of massless colour-sextet quarks is expected to be either a 
Walking or a Conformal field theory. (First term in the j8 function is negative, second positive.) For 
Nf = 3 conformal behaviour is expected. Nf = 2 could, a priori, exhibit either behaviour. We use 
unimproved lattice QCD with staggered quarks to simulate these theories. The RHMC simulation 
algorithm is used to allow us to simulate at values of Nf which are not multiples of 4. 

We simulate Nf = 2 thermodynamics in an attempt to distinguish whether it is walking or 
conformal. Simulations on Nt =4,6,8 lattices show widely separated deconfinement and chiral 
symmetry restoration transitions, which move to weaker couplings as Nt is increased ||8|, ^. For the 
chiral transition, the separation between the N = S and N = 6 transitions is considerably less than 
that between the Nt = 6 and Nt = 4 transitions. 2-loop perturbation theory suggests that between 
Nt = S and N = 6 the shift is consistent with asymptotic freedom, while between N = 6 and N = 4 
the theory is strongly-coupled and the shift is due to quenched evolution of the coupling. If this 
is correct, the deconfinement transitions occur in the strong-coupling (quenched) regime for any 
Nt we could reasonably study. We are thus limiting our studies to the evolution of the couplings 
at the chiral transition as N is increased. Other groups are studying QCD with 2 colour-sextet 
quarks using different lattice actions [ |T0| , pi] , 12, |l^, 14, 15, |16|, |l^]. Most of these studies are of 



the zero-temperature behaviour of this theory. A consensus as to whether this theory is QCD-like 
or conformal has yet to be reached. 

We are now extending our simulations to Nt = 12 lattices and plan simulations on A', = 18 
lattices, to check if this theory is indeed QCD-like, or whether we are observing the approach to a 
bulk chiral transition, expected for a conformal theory. (Our simulations aie limited to the neigh- 
bourhood of the chiral transition - weak coupling.) Since, however, 2-loop perturbation theory 
predicts that Pxi^f ~ ~ Pxi^f = 8) ~ 0.12, we first need to increase the number of j8 values 
and statistics in the neighbourhood of the chiral transition, for Nt = S. We are also adding simula- 
tions at a smaller quark mass, m = 0.0025. These additional N = S runs also provide evidence that 
the Nf = 2 chiral transition is second order. 



2. Simulations at A/^ = 8 

We are extending our simulations on 16^ x 8 lattices in the neighbourhood of the chiral tran- 
sition. Our quark masses are m = 0.0025, m = 0.005, m = 0.01, m = 0.02, which will allow us to 
access the chiral limit. Whereas our earlier simulations used jS values separated by 0.1, we have 
decreased this to 0.02 close to the transition. Near this transition we use 50,000-trajectory runs 
at each (j8,m) for m = 0.0025, 20,000-trajectory runs for m = 0.005 and m = 0.01 and 10,000- 
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trajectory runs at m = 0.02. We are currently increasing or plan to increase our run lengths at each 
of these masses. 

The chiral condensates for each mass decrease as jS increases. More importantly, as jS in- 
creases, the mass dependence of these condensates becomes more pronounced. The decrease in the 
chiral condensate with decreasing mass is such that it does appear that it will vanish in the chiral 
limit for jS sufficiently large. However, the j3 dependence of {^f^lf) is sufficiently smooth at all the 
masses of our simulations, that we would need a precise analytical form to perform a believable 
chiral m -> extrapolation to determine where it vanishes. This we do not have. Hence we examine 
the (disconnected) chiral susceptibilities 



(2.1) 



where V is the spatial volume of the lattice and T is the temperature. is a lattice averaged 
quantity. Because we only have stochastic estimators for (5 per trajectory), we obtain un- 
biased estimators of (yy)^ as the products of 2 different estimators of ij/Y for the same gauge 
configuration. 
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Figure 1: Chiral susceptibilities on a 16^ x 8 lattice. 



The chiral susceptibility diverges at the chiral phase transition for zero quark mass. At small 
but finite mass, it shows a clear peak which becomes sharper as m decreases. Extrapolating the 
position of said peak to m = yields P^, the j8 value of the chiral phase transition. Figure |l| shows 
the chiral susceptibilities from our runs on 16^ x 8 lattices. What is clear from this plot is that the 
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position of the peak in the chiral susceptibility has very little dependence on the quark mass m. It 
is this fact that allows us to extrapolate its position to m = 0. Our best estimate from the 'data' 
presented here is = 6.70(2). 
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Figure 2: Peaks of the chiral susceptibilities for each mass as functions of the scaling variable m'/^^'. 

Assuming that the chiral transition is a second order transition in the equivalence class of the 
3-dimensional 0(2) or 0(4) spin model, as expected for a tinite temperature transition, we fit the 
peaks of the susceptibilities to the scahng form 

Xiif^^{max)=Cm^l^-\ (2.2) 

with 8 = 4.8. Figure ^ shows that the fit is quite reasonable - jDOF 2.2. This contrasts to the 
fit to the form for a first-order transition (5 = °°) which is very poor - "x} jDOF 65. 

3. Simulations at A/f = 12 

Our A^; = 12 simulations are being run on a 24^^ x 12 lattice for quark masses m = 0.01, m = 
0.005 and m = 0.0025 and a range of j8 values which includes the chiral transition. We perform runs 
of 10,000 length- 1 trajectories away from this transition. Close to the transition (6.6 < j8 < 6.9) the 
j8 values are spaced by 0.02, and we plan runs of 50,000 trajectories for each (j3,m). The results 
we are reporting have 10,000-25,000 trajectories at each (j8,m). 

Figure ^ shows preliminary measurements of the chiral susceptibilities obtained in these sim- 
ulations. As one can see, the susceptibility rises rapidly as jS is decreased from j8 = 7.2 down to 
j3 just above 6.8, for each of the 3 mass values. Below this each of the graphs flattens out. This is 
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Figure 3: Chiral susceptibilities on a 24^ x 12 lattice. 



7.2 



obvious for m = 0.01,0.005, and there is some preliminary evidence for this at m = 0.0025. From 
this preliminary 'data' we cannot estimate the position of the chiral transition. For this we will 
need higher statistics and possibly a lower mass value. 

Since we cannot as yet estimate the position of the chiral transition from the chiral suscepti- 
bilities, we perform a direct comparison between the chiral condensates on the 16^ x 8 and 24^ x 12 
lattices for the same mass. Ideally we should select the smallest mass {m = 0.0025). However, we 
do not yet have enough j3 values at m = 0.0025 to make a good comparison. Hence we compare 
the chiral condensates at m = 0.005 on the two lattices. Figure |^ shows this comparison. The 
24^ x 12 condensates are displaced to larger j8 with respect to the 16^^ x 8 condensates. However, 
the displacement appears to be less than the « 0.12 predicted from the 2-loop Callan-Symanzik 
beta-function. 



4. Discussion and conclusions 

We are simulating lattice QCD with 2 flavours of light staggered quarks at finite temperatures 
on 16^ x 8 and 24^ x 12 lattices for j3 close to the chiral transition. Our goal is to determine whether 
the chiral transition for massless fermions approaches zero coupling (infinite j3) as A/^, ^ oo in a 
manner predicted by the perturbative j3 -function, implying a QCD-like (walking) theory, or if it 
approaches a finite non-zero value indicating a conformal field theory. For A/^r = 8 we are able to 
determine the position of the chiral transition with some accuracy (j3;^ = 6.70(2)) from the peaks 
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Figure 4: The m — 0.005 chiral condensates on 16^ x 8 and 24^ x 12 lattices as functions of j3. 



in the chiral susceptibilities. The scaling of these peaks with mass is consistent with a second-order 
phase transition in the 0(2) / 0(4) universality class, expected of a finite-temperature transition, and 
inconsistent with the first-order transition expected for a bulk transition. For Nt = 12 our 'data' is 
still too preliminary to determine from the chiral susceptibilities. However, direct comparison of 
the chiral condensates suggests that ^^{Nt = 12) > j3;(.(A'f = 8) as predicted by asymptotic freedom 
for a finite temperature transition. 

For Nt = 8, we have also performed simulations on a 24^ x 8 lattice at a few strategically 
chosen values of (j8,m). Comparison of the chiral condensates and chiral susceptibilities between 
the 2 spatial volumes indicates that the finite volume effects are small. 

For A/, = 12, we need more statistics at more values of j8 for each mass, to determine the peaks 
in the chiral susceptibilities. If the susceptibilities remain too flat to be able to determine their 
peaks, we may need to simulate at an even smaller mass (m = 0.00125) to determine We will 
also need to check that the finite volume effects are under control. We see no reason why results as 
precise as those we obtained on A/, = 8 lattices cannot be obtained. 

It is probable that simulations at even larger Nt& will be needed to resolve the true nature of 
this theory. We are planning simulations at A/^r = 18 (36^ x 18 lattices). In addition we are planning 
some zero temperature simulations in the chirally-restored phase. 
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